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ABSTRACT 

Recent findings by 7-ray Cherenkov telescopes suggest a higher transparency of the 
Universe to very-high-energy (VHE) photons than expected from current models of 
the Extragalactic Background Light. It has been shown that such transparency can be 
naturally explained by the DARMA scenario, in which the photon mixes with a new, 
very light, axion-like particle predicted by many extensions of the Standard Model 
of elementary particles. We discuss the implications of DARMA for the VHE 7-ray 
spectra of blazars, and show that it successfully accounts for the observed correlation 
between spectral slope and redshift by adopting for far-away sources the same emission 
spectrum characteristic of nearby ones. DARMA also predicts the observed blazar 
spectral index to become asymptotically independent of redshift for far-away sources. 
Our prediction can be tested with the satellite-borne Fermi/LAT detector as well 
as with the ground-based Cherenkov telescopes H.E.S.S., MAGIC, CANGAROOIII, 
VERITAS and the Extensive Air Shower arrays ARGO-YBJ and MILAGRO. 
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, An impressive amount of information about the Uni- 

■ verse at 7-ray energies larger than 100 GeV - namely in 
the very-high-energy (VHE) 7-ray band - has been col- 

,, lected over the past few years by the Imaging Atmospheric 

■ Cherenkov Telescopes (lACTs) such as H.E.S.S., MAGIC, 
CANGAROO III and VERITAS. These lACTs have de- 
tected gamma-ray sources over an extremely wide interval 
of distances, ranging from the parsec scale for Galactic ob- 
jects up to the Gpc scale for the farthest blazar 3C 279 at 
redshift 2 = 0.536 (Albert et al. 2008). 

These observations allow both to infer the intrinsic 
properties of the sources and to probe the nature of pho- 
ton propagation through cosmic distances. The latter fact 
becomes particularly important in connection with VHE 7- 
ray observations, since hard photons travelling through cos- 
mological distances interact with soft background photons 
permeating the Universe, producing e"'"e~ pairs through the 
standard 77 e^e~ process and thereby disappearing. De- 
noting by E and e the energy of the hard (incident) and 
of the soft (background) photon, respectively, and by tp the 
scattering angle, the corresponding cross section is (Heitler 
1960) 

(j^^{E,e,ip) = 1.25 X 10"^^ (1 - /3^) x 

X 2(3{0'-2) + (3-/3*)lni±| cm% (1) 



which depends on E, e and ip through the dimensionless 
parameter* l3{E,e,ip) = [l — 2me/-Ee (1 — cos(^)] This 
process is kinematically allowed for e above the energy 
threshold ethr(£^,V') = 2mg/i? (1 — cos (/p), i.e. for /3 > 0. 
ajj{E,e,(p) reaches its maximum, aijiy^" ~ 1.70 ■ 10~^^cm^, 
for P ~ 0.70. Assuming head-on collisions, a^^^ is at- 
tained when the background photon energy is e*(-E) — 
(0.5TeV/i?) eV. This shows that in the energy interval 
100 GeV < £ < 100 TeV, explored by the lACTs, the re- 
sulting opacity is dominated by the interaction with in- 
frared/optical/ultraviolet diffuse background photons - usu- 
ally called Extragalactic Background Light (EBL) - with 
0.005 eV < e < 5eV (corresponding to the wavelength range 
0.25 /im < A < 250 /im). 

The EBL is produced during the star-formation history 
of the Universe, including a possible early generation of stars 
formed before galaxies were assembled. Based on synthetic 
models of the evolving stellar populations in galaxies as well 
as on deep galaxy counts (see, for a review, Hauser & Dwek 
2001), several estimates of the spectral energy distribution 
(SED) of the EBL have been proposed, leading to different 
values for the transparency of the Universe to 100 GeV < 

* Natural Lorentz-Heaviside units with h = c = = 1 are 
employed throughout. 
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E < 100 TeV photons (Stecker et al. 1992 and 2006; Knciske 
et al. 2002 and 2004; Primack et al. 2005; Aharonian et al. 
2006; Stanev et al. 2006; Mazin & Goebel 2007; Franceschini 
et al. 2008); the resulting uncertainties are large. Besides 
including the evolution arising from the cosmic expansion, 
the evolutionary effects of the stellar populations of galaxies 
should be taken into account in the evaluation of the SED 
of the EBL (Raue & Mazin 2008). 

Because of the absorption produced by the EBL, the 
propagation of a monochromatic photon beam of observed 
energy Eq from a source at redshift z is controlled by the op- 
tical depth T^{Eo, z)^. Therefore, the observed photon spec- 
trum $obB(-Bo,x) is related to the emitted one ^ein{E(z)) 
by 

<i>ob.{Eo, z) = e-^^^^"'^) $em {Eq{1 + z)) . (2) 

Clearly t^{Eo,z) increases with z, since a greater source 
distance entails a larger probability for a beam photon to 
be absorbed. It is generally assumed that $em oc E~^'"° for 
E > 100 GeV, and so from Eq. (2) it follows that 3>obs(-Eo, z) 
decreases exponentially with z. From now on, we restrict 
our attention to the energy range 0.2 TcV < Eq < 2 TeV in 
which the blazar spectra have been measured by lACTs. In 
this range, data axe fitted as $obs oc E'^^^". Hence, Eq. (2) 
yields 

r„bs{z) ^re^ + T^{Eo,z) , (3) 

up to a logarithmic ^-dependence. The experimental results 
are summarized in Table 1 and the spectral slopes are plot- 
ted vs. redshift in Fig.l. 

A full-fledged prediction of how a broadband 
(100 GeV < .Bo < 100 TeV) SED of blazars is affected by the 
EBL would be beyond the scope of this paper. We will, in- 
stead, examine modifications of blazar spectra duo to physi- 
cal processes possibly taking place in intergalactic space and 
discuss them with respect to current lACT data. Such data 
cover the 0.2 TeV < Eq < 2 TeV energy range, hence the 
corresponding relevant EBL range is 0.25 eV < eo < 2.5 eV. 
Within this range, we find it convenient to adopt the follow- 
ing analytic parametrization of the EBL spectral number 
density at the present cosmic epoch (Stecker et al. 1992, 
SDS) 

(, \ -2-55 
cm-^eV"' , (4) 

where a is a suitable constant. Indeed, in the range 
0.25 eV < eo < 2.5 eV, Eq. (4) does reproduce the basic be- 
havior of the SED reported in the most recent phenomeno- 
logical model of the EBL (Franceschini et al. 2008) for 0.5 < 
a < SJ. We remark that in the band 0.25 eV < eo < 2.5 eV 



t We are adopting a standard ACDM cosmological model with 
f2jv ~ 0.7 and Qm — 0.3. Distances are expressed in terms of the 
redshift z and we shall henceforth write E{z) = Eo(l + z) and 
e(z) = eo{l + z), with Eo and eo referring to the present cosmic 

epoch (2 = 0). 

t The SED in the model of Franceschini et al. (2008) is 
slightly convex in the log [eon^(eoiO)] vs. logeo representation 
for 0.25 eV < eo < 2.5 eV, due to the emission bump resulting 
from the integrated emission of the low-mass star population that 
remains close to the main sequence over cosmological times (see 
e.g. Kneiske et al. 2002). The values a = 0.5 and a = 3 bracket a 



Table 1. Blaaars with known redshift and VHE 7-ray flux and 
spectrum. 
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PKS 0548-322 


0.069 


2.80 ± 0.30 


(3.3±0.7] 


xlO" 


12 


PKS 2005-489 


0.071 


4.00 ± 0.40 


(3.3±0.5] 


xlO" 


12 


RGB J0152-I-017 


0.080 


2.95 ± 0.36 


(4.4±1.2] 


xlO" 


12 


PKS 2155-304 


0.116 


3.37 ±0.07 


(2.9±0.2) 


xlO" 


11 


lES 1426-1-428 


0.129 


3.55 ± 0.46 


(2.5±0.4^ 


xlO" 


11 


lES 0229-1-200 


0.139 


2.50 ± 0.19 


(4.5±0.7; 


xlO" 


12 


H 2356-309 


0.165 


3.09 ± 0.24 


(2.5±0.7; 


xlO" 


12 


lES 1218-1-304 


0.182 


3.00 ± 0.40 


(1.0±0.35 


xlO" 


11 


lES 1101-232 


0.186 


2.94 ±0.20 


(4.4±0.75 


xlO" 


12 


lES 0347-121 


0.188 


3.10 ±0.23 


(3.9±0.7) 


xlO" 


12 


lES lOlH-496 


0.212 


4.00 ± 0.50 


(6.4±0.3) 


xlO" 


12 


PG 1553-H13 


>0.25 


4.20 ± 0.30 


(5.2±0.9) 


xlO" 


12 


30 279 


0.536 


4.1 ±0.7 


(2.9±0.5) 


xlO" 


11 



The observed photon spectral index Fobs in the 0.2 TcV < i?o < 
2 TeV band (0.2 TeV < £0 < 0.6 TeV for 3C 279: the effective 
0.2 TeV < -Eq < 2 TeV slope of this source might be steeper) and 
the flux <I>obs(> 0.2 TeV) measured at energy Eo > 0.2 TeV (in 
ergcm"^ s"'^), for blaaars observed at different redshift z. The 
errors indicate the statistical uncertainty; the corresponding sys- 
tematic uncertainties on the spectral index are typically ~0.1 for 
H.E.S.S. and ~0.2 for MAGIC. See De AngeUs, Mansutti & Persic 
(2008) for references and more informations. 

the assumed SDS parametrization also encompasses both 

the lowest EBL level predicted by Primack et al. (2005) and 
the higher EBL level predicted by the "baseline model" of 
Stecker et al. (2006)§. 

From the definition of r-^(So, z) (see e.g. Fazio & Stecker 
1970), we finally obtain 

r,(i?o,.)=. 0.50a (^^^)"'' [(1 + ^)*" " l] , (5) 

where evolutionary effects arising from galaxy evolution 
(Raue & Mazin 2008) have been included on top of those 
produced by the cosmic expansion. Wo have chocked that 
for 1.2 < a < 3, 0.2 eV < eo < 2.5 cV and 0.05 < z < 0.5 
Eq. (5) is approximately consistent with the analytic fit, pro- 
vided by Stecker & Scully (2006), of T^(i5o,2) as predicted 
by the 'baseline model' of Stecker et al. (2006). 

Since for the nearby blazars (z < 0.03) the EBL photon 
absorption is expected to be negligible, we assume that for 
such sources observations do yield 3>em- We find, on aver- 
age, Fem — 2.4 \ Various emission models for blazars have 

linear stripe in the log [eo ™7(eo, 0)]-log£o plane which envelopes 
the actual EBL shape. 

§ In this paper we disregard the "fast evolution model" of Stecker 
et al. (2006), which predicts an even higher EBL level and would 
require a > 3. 

^ Several interpretations for the scattering of the data in Fig. 1 
for 2 < 0.2 are possible: e.g., different sources can be observed in 
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Figure 1. Observed values of the spectral index for all blazars 
detected so far in the VHE band are represented by dots and 
corresponding error bars. Superimposed on them is the predicted 
behaviour of the observed spectral index within two different sce- 
narios. In the first scenario (light grey area) T^^^^ is computed in 
terms of standard physics in the SDS model of the EBL. In the 
DARMA scenario (dark grey area) rj^j^^*^^^ is evaluated within 
the proposed photon-ALP oscillation mechanism as based on the 
same SDS model of the EBL. 

been developed and are briefly summarized in the Appendix. 
In the widely used Synchro-Self-Compton (SSC) emission 
model (e.g., Ghisellini et al. 1998) Tom — 2.4 suggests a 
Compton peak at around or below 100 GeV. The predicted 
observed spectral index r^^g^ then follows from Eqs. (3) and 
(5): it is represented as the light-grey area in Fig. 1. 

Fig. 1 shows that the actually observed spectral index 
increases more slowly than V^^^^ for redshifts z > 0.2. More- 
over, the observed values cannot be explained for z > 0.3 
by the EBL model of SDS even for a as low as 0.5. Being 
the optical depth a monotonically increasing function of z, 
we interpret the conflict between Fobs and F^^g^ - apparent 
in Fig. 1 - as calling for a departure from the conventional 
view and we proceed to explore its consequences. 

A possible way out relies upon a systematic harden- 
ing of the emission spectrum with increasing z. That is. 
Fern - which is currently supposed to be independent of z - 
has to decrease as z increases, so as to offset the growth of 
T.y{Eo,z). This situation is very difficult to achieve within 
the SSC emission model: even assuming that we selectively 
observe increasingly flaring sources at higher redshifts, the 
radiating electrons will be emitting more and more in the 
Klein-Nishina regime. Therefore - unlike the synchrotron 
peak which will appreciably shift to higher energies - the 
Compton peak will hardly shift, thereby ensuring that Fem 
is indeed virtually independent of z (Persic & De Angelis 
2008). 

Other possibilities have been suggested, based on mod- 
ifications of the emission mechanism. One involves strong 
relativistic shocks, which can give rise to values of Fem con- 
siderably smaller than previously thought (Stecker et al. 
2007; Stecker & Scully 2008). Another rests upon photon 
absorption inside the blazar, which has be shown to lead 

different emission states, and so they may exhibit slightly different 
slopes (Persic & De Angelis 2008). 



to a substantial change of Fem (Aharonian et al. 2008a). 
A further option could be the underlying emission mecha- 
nism of the more luminous sources (e.g., flat-spectrum radio 
quasars) being hadronic, with the muon (and cascade) syn- 
chrotron component peaking in the sub-TeV region and the 
(neutral- and charged-)pion cascades crossing the 0.2-2 TeV 
band with a very hard (F ~ 1.9) spectrum, as shown in the 
very case of 3C 279 (Bottcher et al. 2008). 

While successful at substantially reducing Fem in indi- 
vidual sources, these attempts fail to provide a systematic 
in the correlation of Fem versus z, needed to overcome the 
above difficulty. 

In the framework of the same basic emission mechanism 
being at work in blazars and flat-spectrum quasars both at 
low and at high redshift, we propose a different solution, re- 
ferred to as the DARMA scenario (De Angelis, Roncadelli 
& Mansutti 2007). Implicit in all previous considerations is 
the hypothesis that photons propagate in the standard way 
throughout cosmic distances. We suppose instead that, in 
the presence of cosmic magnetic fields, photons can oscillate 
into a new, very light, spin-zero particle - named Axion-Like 
Particle (ALP) - and vice-versa. Once ALPs are produced 
close enough to the source, they travel unimpeded through- 
out the Universe and can convert back to photons before 
reaching the Earth. Since ALPs do not undergo EBL absorp- 
tion, they act as if the observed photons had an effective op- 
tical depth smaller than Tj{Eo,z). Eq. (3) entails that Fobs 
gets reduced by the same amount as far as the ^-dependence 
is concerned, thereby avoiding the conflict shown in Fig. 1. 
Now the dependence of Fobs on z is in agreement with ob- 
servations, since the photon-ALP oscillation reduces photon 
absorption even for standard emission spectra. In order to 
guarantee consistency with observations of nearby blazars 
we take Fem — 2.4 for all sources represented in Fig. 1. 

The key ingredient of the DARMA scenario - namely 
the existence of ALPs - is not an ad hoc assumption in- 
vented to solve the problem in question|l. Instead, very light 
ALPs turn out to be a generic prediction of many extensions 
of the Standard Model of elementary particle physics and 
have attracted considerable interest over the past few years. 
Besides than in four-dimensional models (Masso & Toldra 
1995 and 1997; Coriano & Irges 2007; Coriano et al. 2007), 
they naturally arise in the context of compactified Kaluza- 
Klein theories (Chang et al. 2000; Dienes et al. 2000) as well 
as in superstring theories (Turok 1996; Svrcek & Witten 
2006). Moreover, it has been argued that an ALP with mass 
10"^^ eV is a good candidate for the quintessential 
dark energy (Carroll 1998) which might trigger the present 
accelerated cosmic expansion. 

Below, we first review the motivation and the properties 
of ALPs which are of direct relevance for the present discus- 
sion. Next, we outline the computation of the predicted ob- 
served spectral indices FJ^,^^'^^. Details on the derivation, 
as well as on the dependence of FJ^j^^^^ on the adopted 
EBL model, will be reported elsewhere. 



II Other aspects concerning the relevance of ALPs for gamma- 
ray astrophysics have been addressed in (Hooper & Serpico 2007; 
Hochmuth & Sigl 2007; De Angelis, Mansutti & Roncadelli 2008; 
Simet et al. 2008; Cheleuche et al. 2008, Cheleuche & Guendelman 
2008). 
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On the basis of phenomenological as well as concep- 
tual arguments, the Standard Model is currently viewed as 
the low-energy manifestation of some more fundamental and 
richer theory of all elementary-particle interactions includ- 
ing gravity. Therefore, the lagrangian of the Standard Model 
is expected to be modified by small terms describing inter- 
actions among known and new particles. ALPs are spin-zero 
light bosons defined by the following low-energy eflfective lar- 
grangian 

jCal.p = ^d*" ad^ia-^m^ - F'"" a , (6) 

where F^" is the electromagnetic field strength, F^,^ is its 
dual, a denotes the ALP field whereas m stands for the ALP 
mass. According to the above view, it is assumed for the 
inverse two-photon coupling M » Gp ~ 250 GeV. On 
the other hand, it is supposed that m ^ Gp ~ 250 GeV 
and for dcfinitencss wc take m < 1 eV. As far as generic 
ALPs are concerned, the parameters M and m are to be 
regarded as independent. 

So, what really characterizes ALPs is the trilinear 7-7-a 
vertex described by the last term in £alp, whereby one ALP 
couples to two photons. Owing to this vertex, ALPs can be 
emitted by astronomical objects of various kinds, and the 
present situation can be summarized as follows. The nega- 
tive result of the CAST experiment designed to detect ALPs 
emitted by the Sun yields the bound M > 0.86-10"' GeV for 
m < 0.02 eV (Zioutas et al. 2005; Andriamoje et al. 2007). 
Moreover, theoretical considerations concerning star cooling 
via ALP emission provide the generic bound M > GeV, 
which for m < 10"^'^ eV gets replaced by the stronger one 
M > 10" GeV even if with a large uncertainty (Raffelt 1990 
and 1996; Khlopov et al. 2004). 

The same 7-7-a vertex produces an ofi-diagonal element 
in the mass matrix for the photon- ALP system in the pres- 
ence of an external magnetic field B. Therefore, the interac- 
tion eigenstates differ from the propagation cigenstates and 
photon- ALP oscillations show up (Sikivie 1984a and 1984b; 
Maiani et al. 1986; Raffelt & Stodolsky 1988). The situation 
is analogous to what happens in the case of massive neu- 
trinos with different fiavours, but while all neutrinos have 
equal spin - hence neutrino oscillations can freely occur - 
ALPs have instead spin zero whereas the photon has spin 
one, and so the transformation can take place only if the 
spin mismatch is compensated for by an external magnetic 
field B. 

We imagine that a sizeable fraction of photons emitted 
by a blazar soon convert into ALPs. They propagate un- 
affected by the EBL and we suppose that before reaching 
the Earth a substantial fraction of ALPs convert back into 
photons. We further assume that this photon- ALP oscilla- 
tion process is triggered by cosmic magnetic fields, whose 
existence has been demonstrated very recently by AUGER 
observations (Abraham et al. 2007). Lacking information 
about their morphology, we assume that cosmic magnetic 
fields have a domain-like structure (Kronberg 1994; Grasso 
& Rubinstein 2001; Furlanetto & Loeb 2001): B is supposed 
to be constant over a domain of size Ldom equal to its coher- 
ence length and to change randomly its direction from one 
domain to another while keeping approximately the same 
strength. As argued by De Angelis, Persic & Roncadelli 
(2008), it looks plausible to assume the coherence length 



in the range 1 Mpc < Ldom < 10 Mpc. Correspondingly, the 
inferred strength lies in the range 0.3 nC < Bo < 0.9 nG. 
This conclusion agrees with previous upper bounds (Blasi 
et al. 1999), so we assume Ldom = 7 Mpc and Bo = 0.5 nG 
as reference values at z = 0. 

Following the same computational procedure as in De 
Angelis, Roncadelli & Mansutti (2007), we evaluate the 
probability Py^^{Eo,z) that a photon remains a photon 
after propagation from the source to us when allowance is 
made for photon- ALP oscillations as well as for photon ab- 
sorption by the EBL. As a consequence, Eq. (2) becomes 

$obs(So, Z) = P-y^-yiEo, Z) $em (Fo(l + z)) (7) 

so that Eq. (3) gets replaced by 

robs(2) ~rem-lnP^^^(So,2) , (8) 

again up to a logarithmic z-dependence. Assuming m < 
10~^° eV and M ~ 4- lO" GeV as in De Angelis, RoncadeUi 
& Mansutti (2007) and adopting the same EBL model as be- 
fore (see Eq. (4) with 0.5 < a <3), our result for Pot'^'^'^^ is 
represented by the dark gray area in Fig. 1. Wc have checked 
that the same result remains practically unaffected within 
the range 10^^ GeV < M < 10^^ GeV. 

In conclusion, for a realistic EBL model (defined by 
Eq. (4) with 0.5 < a < 3) and assuming the same nomi- 
nal emission spectral slope Pom — 2.4 for all VHE blazars, 
the DARMA scenario naturally explains the lACT data 
and predicts that Fobs becomes asymptotically indepen- 
dent of z for far-away sources. Our prediction can be tested 
with the satellite-borne Fermi/LAT detector as well as with 
the ground-based Cherenkov telescopes H.E.S.S., MAGIC, 
CANGAROO III, VERITAS and the Extensive Air Shower 
arrays ARGO-YBJ and MILAGRO. We remark that the 
DARMA scenario could loose much of its motivation - and 
be eventually disproved - if the emission mechanisms of 
VHE blazars and quasars had a variation according to, 
e.g., luminosity. The most distant and luminous VHE 7- 
ray source that appears in Fig.l is 3C279, a remarkable 
example of fiat-spectrum radio galaxies: for these sources, 
the flaring and accompanying intermittency of source activ- 
ity, not well understood at present, may point to emission 
mechanisms different from those that are commonly being 
used for blazars (i.e., the leptonic SO models). Such emis- 
sion mechanisms may provide the hard spectra emitted by 
high-luminosity, high- 2: sources, which in Fig.l are required 
to counterbalance the spectral steepening imposed to TeV 
radiation by traversing the EBL over cosmological distances. 
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APPENDIX 

The observed emission from a blazar is mainly interpreted 
as radiation from a relativistic jet that is directed approx- 
imately along the line of sight to the observer (Urry & 
Padovani 1995). 
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Leptonic models associate the TeV emission with the 
synchrotron-Compton (SC) process occurring in the jet. In 
one/more magnetized plasma blob(s) moving relativistically 
toward the observer, the electrons emit synchrotron radia- 
tion that is, in turn, Compton upscattcrcd by their par- 
ent electrons (one/multi-zone synchro-self-compton [SSC] 
scheme: e.g., Maraschi et al. 1992). A variant of this scheme 
involves an additional, external field of soft photons - e.g., 
the infrared field associated with the broad-line region (ex- 
ternal inverse Compton [EIC] scheme: Dermer ct al. 1992). 
Such leptonic models are mainly motivated by the double- 
hump shape generally observed in blazar SEDs. However, 
looscncss/lack of X-ray/TcV variability correlation includ- 
ing 'orphan' TeV flares (e.g., Krawczynski ct al. 2004; Blaze- 
jowski ct al. 2006), and peculiar quasar SEDs (e.g., Bottcher 
et al. 2008) increasingly challenge the simple - hence popu- 
lar - SC model. 

Hadronic processes arc more complicated: they can 
entail photon-initiated cascades (Mannheim & Biermann 
1992) or pp collisions (Pohl & Schlickeiser 2000; Dar & Laor 
1997; Beall & Bednarek 1999), but both essentially produce 
VHE photons via tt" decay. For TeV blazars it was suggested 
that synchrotron emission from ultra-relativistic protons in 
strong magnetic fields is mostly responsible for the observed 
TeV emission from blazars (Aharonian 2000; Miicke et al. 

2003) . The hadronic model may successfully describe the 
broad-band SED ofblazars (Aharonian et al. 2005; Bottcher 
et al. 2008) and, in principle, account for the X-ray/TcV cor- 
relation, if an appreciable amount of X-ray emission comes 
from the synchrotron radiation of secondary electrons. There 
is in fact some flexibility here, as the co-accelerated, primary 
electrons might also contribute to the X-ray band in a sig- 
nificant manner. This possibly 'double' origin of the X-rays 
may explain the scatter in the X-ray/TeV correlation. 

Current leptonic and hadronic models, however, are 
challenged by the 'orphan' flares and by the minute/hour- 
scalc variability of blazars both at TeV energies (Gaidos et 
al. 1996) and X-ray energies (Cui 2004; Xue & Cui 2005). 
The flaring and accompanying intermittency of source ac- 
tivity is probably a key to understanding the physics of flat- 
spectrum radio sources, of which the outlier in Fig.l, 3C 279, 
is a remarkable example. Furthermore, if ultra-high-energy 
{E> 10^" eV) cosmic rays (Linsley 1963; Protheroe & Clay 

2004) are produced in AGNs, then purely leptonic models 
are certainly limited and - eventually - inadequate. For gen- 
eral reviews on the origin of astropliysical TeV photons, see 
e.g. Cui (2006), Aharonian et al. (2008b), and DeAngelis, 
Mansutti & Persic (2008), and references therein. 
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